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Abstract -Recent data from Reticulum II (Retll) require the energy range of the FermiLAT 
7 -excess to be ~ 2 — 10 GeV. We adjust our unified nonthermal Dark Matter (DM) model 
to accommodate this. We have two extra scalars beyond the Standard Model to also explain 
3.55 keV X-ray line. Now the mass of the heavier of them has to be increased to lie around 250 
GeV, while that of the lighter one remains at 7.1 keV. This requires a new seed mechanism 
for the 7 -excess and new Boltzmann equations for the generation of the DM relic density. 
All concerned data for Retll and the X-ray line can now be fitted well and consistency with 
other indirect limits attained. 


The endeavour for the detection of Dark Matter (DM) 
is increasingly gaining momentum. Gamma-ray signals 
from the FermiLAT experiment have attracted much 
attention [ m- mi These cannot be explained by the 
known astrophysical processes. On the other hand, their 
DM origin has been a topic of debate [ [T2] ~ [42]]. One 
possibility is the decay/self-annihilation of DM particles 
clustered around massive gravitating bodies, e.g. the 
Galactic Centre (GC) or dwarf galaxies. Separately, 
an X-ray line of energy 3.55 keV has been reported [ 
[43] . [44] ] by the XMM Newton observatory by use of a 
data set obtained from Andromeda and 73 other galaxy 
clusters including Perseus. An astrophysical explana- 

(a) present address: Harish-Chandra Research Institute, Chhat- 
nag Road, Jhusi, Allahabad, 211019, India 


tion [45] of this line, though possible, is beset :46{ with 
uncertainties in the potassium abundance in the target. 
Thus a DM origin of the X-ray line remains a viable 
possibility and could be from decaying [ [T7], [3B3] an¬ 
nihilating m or excited DM. It would be a worth¬ 
while effort to construct a unified DM model for these 
two phenomena. 

Data from the dwarf spheroidal galaxy Retll mi sug¬ 
gest an upward shift in the earlier claimed i®~ 0 ] en¬ 
ergy range of the FermiLAT 7-excess to 2 —10 GeV. The 
high galactic latitude of Retll makes its 7-emission rel¬ 
atively free from complicated backgrounds. This higher 
range is what we adopt here. That requires a modifi¬ 
cation in our 2-component nonthermal DM model m , 
proposed earlier to explain both the 7-excess and the 
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X-ray line. In our model the fields describing DM have 
tiny couplings with Standard Model (SM) fields. As a 
result, the DM particles are produced nonthermally and 
they are unable to thermalise later. Two extra elec- 
troweak (EW) singlet scalar fields S 2 3 are introduced. 
These and the SU(2 )l doublet Higgs field H comprise 
the scalar sector. Inter-mixing among them leads to 
three physical particles x t 2 3 with M x ^ ~ 125 GeV, x 2 
and x (with M x ~ 7 keV) having tiny mixing angles 
between them. The decays X 3 = y 77 and X 2 _■ y bb (with 
the b’s emitting neutral pions via hadronisation) respec¬ 
tively account for the X-ray line and 7 -excess. Relic 
DM, a mixture of X 2 and X 3 > forms after EW sym¬ 
metry breaking through the processes x 
W + W~ X 2 3 X 2 ' 3 , ZZ -x X. „x. „, tt, 

XjXi 


•^■2, 3 ^2, 3 ’ 


X"2, 3 ^2, 3 5 
A-2, 3 ^2, 3 * 


-^2, 3 ^2, 3' 


An important feature here is the sensitive link be¬ 
tween M x and the energy spectrum of the 7 -excess. 
Indeed, we need M x in the ballpark of 250 GeV to 
fit the increased energy range of this excess. As shown 
numerically later, too small a magnitude of M x , as 
compared with this ballpark value, would unacceptably 
shift the energy spectrum of the 7 -excess to a lower 
range. On the other hand, too large a mass of x 2 
would inhibit its pair production which took place after 
the EW phase transition (Tw ~ 153 GeV [5T]). Now 
the decay x, —► X 2 X 2 is disallowed and X 2 ’ s are pro- 
duced in the early Universe from the pair annihilation 
of SM fermions and gauge bosons. Moreover, the de- 
cay X 2 —>■ W + W~ is now allowed. The strength of the 
X 2 W + W~ (x 2 bb) coupling is proportional to M^/v 
(grrib/Mw), g being the SU(2 )l gauge coupling. Con¬ 
sequently, the v —x W + W~ decay channel becomes 
the dominant contribution to our seed mechanism for 
the 7 -excess. 


Let us recount the salient features of our model. The 
stability of all scalar fields is ensured by the discrete 
symmetry Z 2 x Z 2 . With respect to this, S 2 and S 3 
have charges (- 1 , 1 ) and ( 1 , - 1 ) respectively, while those 
of all other SM fields are (1, 1). The scalar potential for 
the Higgs portal is given by V = Vq + V' where 


Vo(H, S 2 ,S 3 ) = Ki (h'H-±v 2 ) 2 + ±k 2 S* 

{S 2 3 -u 2 ) 2 + \plSl 
+\ 12 (H'H)S 2 + A23 s 2 s 2 
+ \ 13 (H'H-lv 2 )(S 2 -u 2 ) , 

V\s 2 ,s 3 ) = aS 2 S 3 . ( 1 ) 


Terms such as (H^H)S 2 S 3 , S 2 S :l . S 2 S 3 are excluded 
by the assumed symmetry. The “small” V' term softly 
and explicitly breaks the Z 2 x Jj ' 2 invariance down to 
that of Z 2 under which S 2 3 are odd and the rest are 
even. This Z 2 is spontaneously broken by the VEV u 
(2MeV < u < lOMevQof S 3 . In © v = (ReH°), H° 
being the neutral member of the doublet 7L, while the 
coupling constants ^ 1 , 2 , 3 , P 2 , Ai 2 , A 2 3 and A 13 obey cer¬ 
tain stability conditions detailed in Ref. m- Domain 
wall formation from the restoration of Z 2 at a high tem¬ 
perature can also be shown to be inconsequential m- 
The physical scalar fields are si = v^Re H° —v, s 2 = 
S 2 and s 3 = S 3 — u with their squared mass matrix 


M 2 


( 2 k 1 v 2 0 2A13 uv 

0 p 2 + \\2V 2 + 2A 2 3 u 2 a 

y 2A13 uv a 2 k 3 u 2 


( 2 ) 


The eigenvalues of m are M 2 with respective eigen- 

*1, 2 , 3 

state fields y . The latter are linearly related to 

^ 1 , 2,3 J 

Sj 2 3 via the mixing angles 6 1 2 , 0 2 3 and $ 13 . These an¬ 
gles are quite tiny because two of them come from sym¬ 
metry breaking also owing to the smallness of the A’s. 
623 is a pure Z 2 x Z 2 symmetry breaking parameter con¬ 
trolled by a which is chosen to be ~ (10 -9 —10 -8 ) GeV 2 
while O 12 is generated by an interplay of a and A 13 which 
has been taken ~ 1CV 9 . The last mixing angle $13 arises 
from the spontaneous breakdown of the Z 2 symmetry 
driven by A 13 . From a UV perspective the smallness of a 
and the A’s could be due to a presumed hidden tree level 
symmetry broken by radiative loop corrections. We fur¬ 
ther choose M x = 125.5 GeV M x ~ 250 GeV and 
A 1 a 2 

My =7.1 keV. The last choice is consistent with an O 

*3 

(MeV) u provided 2 x 10 ~ 7 < K 3 < 4 x 10 -6 . However, 
1 See right panel of Fig island the related discussion. 
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the upper bound is further restricted to ~ 7 x 10“ ‘ if we 
take the X 3 self-interaction cross section a x divided by 
its mass M x to be less than 0.47 cm 2 / g from collisions 
between different galaxy clusters [52], cf. Figd The 
corresponding quantity for x 2 is too small to make any 
difference. 




The decay y i —> X 3 X 3 a ls° contributed. Let Y x be 
the comoving number density (= actual number den¬ 
sity 4- entropy density of the Universe) of y.. It is 

M x 

given as a function of z = —by a set of two coupled 
Boltzmann equations. The latter involve the thermally 
averaged decay width (T Xi x ) as well as the pair- 
production cross section times the relative velocity of 
collision {av) xS ^ for x = W, Z, / and x i>2 - 

Details appear in Ref. m and will not be repeated 
here. The only change is that the decay y -4 X X 2 

dY X 

is disallowed now. Thus, while the equation for d J is 

dY X 

unchanged, that for . 2 is changed to 


dz 


47r 2 


45 x 1.66 


M p i M Xi \Jg+(T) z 2 x 


(EOT, -K." 2 )(ov)„ 

\ CL 


X„ “I" ( CTV )x o X,^X 


2 2 3 3 


(3) 


with a = W, Z, f and y . Further, the DM relic density 
is given (for j = 2, 3) by 

n x h 2 = 2.755 x 10 8 (M x /GeV) Y x (z 0 ), (4) 

where Zo = /To, Tq being the present temperature 
of the Universe. 



10 100 


z(M Z| /T) 

( a ) ^x 2 > I2 x 3 



Figure 2: Feynman diagrams for dominant production chan Figure 3: Variation of relic densities of both the dark matter 
nels of both the dark matter components \ 2 an d X 3 ■ candidates with 2 


Both x 2 an d X got produced nonthermally in the 
early Universe but only after EW symmetry breaking. 
Thereafter, the self-annihilation of W, Z, y and t, (see 
Fig[2]) acted as primary sources of DM particles x 2 3 ■ 


We take as a boundary condition the vanishing of 
Y x at the EW phase transition (z ~ 0.83). Figures 
Ha) and Elb) show the variation of the relic densities 
of both DM candidates v with 2 for different values 

/v 2, 3 
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of A 12 , A 13 which are ~ 10 -9 — 10~ n . Such strengths 
are needed to keep Y x ^ small enough to generate the 
right DM relic density (12 dm h 2 ) at the present epoch. 
Appropriate values have been chosen for A 12 , A 13 de¬ 
pending on whether X 2 or X 3 is the dominant DM com¬ 
ponent. Starting with null values, Y x are seen to rise 
as more and more DM is produced from the decay/self- 
annihilation of SM particles. They eventually saturate 
to respective particular values at z ~ 10 corresponding 
to a temperature T ~ 12 GeV of the Universe, depend¬ 
ing on the particular values of A 12 , A 13 . These satura¬ 
tion values together need to satisfy the PLANCK m 
68 % c.l. constraint 0.1172 < Udm ^ 2 < 0.1226. Con¬ 
tributions from individual pair production channels of 
X 2 towards Q x h 2 are graphically shown in Fig0] with 
chosen parameters given in its legend: blue line for x 1 
green line for Z , red line for W and brown line for t- 
quarks, the last being somewhat less in magnitude. The 
total relic density of X 2 (yellow line) saturates around 
0.06 which is half the total DM relic density (fl-iA 2 ) of 
today, cf. Ref. m- The remainder comes from x 3 ■ 



Figure 4: Contributions of different production channels to 
the relic density of a 300 GeV \ 2 ■ 

The allowed ranges of A 12 , A 13 , A 23 , # 12 , # 13 , ^23 are 
given in Table [T] Given the chosen values of M x and 
u, the range of A 23 is fixed by the need to avoid a late 
time decay of x 2 via x 2 X 3 X 3 ■ The tiny magnitudes 
of 0 i 3 , 023 and 0 i 2 are required by the constraint of keep¬ 
ing the off diagonal elements of A4 2 in @ to be very 
small. Further, the couplings of X 2 with x , which are 
functions of the three A’s and the three 0’s m , remain 


Al2 

^13 

A 23 

~ (1.4-4.5) 

~ (0.7- 2.4) 

~ (0.2-4.7) 

xicr 11 

xlO -9 

xlO -6 

012 

013 

023 

(rad) 

(rad) 

(rad) 

~ (0.001 - 1.75) 

~ (0.063-6.3) 

~ (0.08-1.31) 

xlO" 25 

xl0~ 12 

xlO" 13 


Table 1: Allowed ranges of concerned couplings and mixing 
angles. Also, a ~ (10~ 9 — 10 -8 ) GeV 2 . 


sufficiently feeble to keep the former beyond the reach of 
DM direct detection experiments [ [54]- [55]] . Another 
point to note is that X 2 behaves like a feebly interact¬ 
ing massive particle (FIMP) starting with a vanishing 
number density. Its fractional relic density saturates af¬ 
ter increasing initially (cf. Fig]3jt) as the temperature 
falls in the cooling Universe. This is the hallmark of 
a “freeze-in” behaviour [5B], as contrasted with that of 
a WIMP; the relic density of the latter starts from an 
equilibrium nonzero value, decreases and then freezes 
out at a saturation level. Though much lighter, x 3 also 
freezes in a way similar to that of x 2 (cf. FiglSJa). 

We turn next to the 7 -excess observed from Retll cov¬ 
ering the range 2 — 10 GeV of the FermiLAT 7 -energy 
spectrum. With M x ~ 250 GeV, x 2 — on account 
of its nonzero mixing with the SM-like Higgs boson x 
— decays predominately into W + W~. Because of the 
small x i x 2 X 2 coupling, x 2 pair-annihilation into the 
same final state, via s-channel x x exchange, is a negli¬ 
gible competitor. Ours is the first model explaining the 
Retll 7 -excess from the decay x 2 —► W + W~ with 7 - 
rays coming predominantly out of neutral pions hadro- 
nising from decaying into qq' pairs. Consider the 
7 -flux from Retll at a line of sight distance s and sub¬ 
tending a solid angle All. The differential distribution 
is 


g?<1> 7 

dYldE 


1 JT , 

4- M. dE 

A-2 


(5) 


dN w 

Here , Z is the energy distribution of each of the two 
7 ’s of energy E produced from the W pair, taken nu- 
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merically from Ref. J = 272 represents an average 
of the “astrophysical factor” J [55] over the opening 
solid angle AD = 27t(1 — cosai nt ), the integration angle 
dint being 0.5° [12]. Further, 


J = 


p(s , D)ds dfl , 


( 6 ) 


where p( s, D) describes the variation of the local dark 
matter density in the neighbourhood of Retll. J has 
been taken to be 10 188 GeV cm -2 from Ref. [55] . 
Finally, T' _>.w+w- is the product of the partial 
width for the decay \ 2 W + W~ and the fractional 
relic density for the component x 2 > be. r' x ^.ww = 

-TT7 V x 2 ^ww- Its occurrence in ([5]) is necessitated 
by the two component nature of our DM. The partial 
width, mentioned above, is given in a transparent nota¬ 
tion by 


o 2 

ilWWx a , o 0,1/0 

r X2 ^ +w _ = -^m 8 2 (i - diw^M- 2 ) 1 / 2 x 
M^( 1 - + 12M^M" 4 ) (7) 

with the coupling g ww given by 
2 


2 M 2 / 
v 


(sin 012 cos 023 + COS 012 sin 0 2 3 sin 0 i3 ) 


i — — 

with v = 2 ~*G f 2 , Gf being the Fermi constant. 

The 7 -flux, computed from ([5]), © and 0 for each 

of the three different values of M x , is plotted in FigJS] 

in comparison with the data points. The background 

7 -flux [12] (turquoise line) is also shown. Though the 

computed plots have been generated with T' w+w _ 

x 2 

fixed at 6.27 x 10 -27 s _1 , the fit does not change much, 
as seen by varying the latter through ±0.94 x 10 -27 s _1 . 
In order to produced the above mentioned range of val¬ 
ues of r' the soft Z 2 x Z ' 2 symmetry breaking 

x 2 

parameter a needs to be in the range 10 ~ 9 GeV 2 ^ 
a ^ 10 -8 GeV 2 . Clearly, the fit is worse when M x be¬ 
comes 200 GeV. We have not extended our fits to cover 
X 2 much beyond 300 GeV since the production of x 2 
(say from tt at the EW transition temperature ~ 153 
GeV) is then cut off by phase space. 

Let us discuss indirect constraints on T' w+w _ 
from other observations. First, consider the limit from 



Background + decay of X 2 when 200 GeV- 

Background + decay of x 2 when M fe = 250 GeV . 

Background + decay of x 2 when M fe =300GeV - 

Background 

y-ray Flux observed from Reticulam II ■—>—1 


Figure 5: Energy distribution of the signal for three different 

M x ’s. 

X 2 


the positron flux in the AMS-02 data m- Using this 
data and assuming a single component DM, Ibarra et 
al. [60] plotted a lower limit (their Fig.3) on the partial 
lifetime T~^ J[ ^ W+W _ of the DM particle decaying into 
W + W~ as a function of the DM mass. Since we have 
a two-component DM in our scenario, we need to con¬ 
sider T' x _^ w+w _ instead of F x ^^ w + w ~. (Note that 
the latter reduces to the former when fl x /Dx = 1 i.e. 
one has a single component DM scenario.) We con¬ 
vert the results of Ref. m into a plot of the upper 
limit on _^ w+w _ as a function of the x 2 fractional 
relic density D^/Dx for M x ^ = 250 GeV, 300 GeV. 
These plots are shown in Fig[ 6 ] Note that our cho¬ 
sen value of 6.27 x 10~ 27 s _1 for F^ _^. w + w -, made in 
order to fit the data from Retll, is below (cf. Fig0 
the range of this upper bound so long as D x /Dx is 
less than ~ 0.65(0.9) for M^ = 250GeV(300 GeV). 
Therefore, D x in our model is constrained to be less 
than ~ 0.65 (0.9) times the total relic density Six for 
M %2 = 250 GeV (300 GeV). 

We next turn to the ANTARES [0i null result on the 
observation of muon neutrinos and antineutrinos from 
DM processes at the Galactic Centre. They derived a 
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Figure 6 : r' x _^ w + w - plotted against /Qt’- AMS-02 
upper bound (red for M x = 300 GeV and green for M x ^ = 
250 GeV) as well as our fixed value (black line). 

90% c.l. upper bound on the total flux <f> IV+P(i as a 
function of the mass of the DM particle taking its pair- 
annihilation into bb as the dominant subprocess. This 
is reproduced in the left panel of Fig[7] In our case the 
dominant subprocess is the decay X 2 —1► W + W~. The 
muon neutrino plus antineutrino flux from Retll, con¬ 
sequent upon the decays of the IP’s, is plotted against 
M x in the right panel of Fig[7] Evidently our flux, be¬ 
ing several orders of magnitude lower, is well within the 
ANTARES limit. 




Figure 7: (a) ANTARES upper bound on (from DM 

pair-annihilation), (b) from \ 2 —> W + W~ in our 

model. 

The 3.55 keV X-ray line comes from one of the two 
monoenergetic photons into which x :l decays through 
its tiny mixing with the SM-like Higgs boson ■ The 
corresponding modified partial decay width T x _>, 77 = 


r X;j is constrained to be in the range 2.5 x 10 ~ 29 

s -1 — 2.5 x 10 -28 s ' 1 in order to fit the observed 
data. The computation of r x3 _> 77 is detailed in Ref. 
r‘T?T| and need not to be repeated here. The left (right) 
panel of FigJ5] shows the region in the u — A 13 (u — a) 
plane allowed by the observational constraints. The red 
coloured patch in the left panel is the region compati¬ 
ble with observed 7 -ray and X-ray fluxes as well as the 
PLANCK limit on the total DM relic density. Similar 
is the case with the patch in the right panel. It is clear 
from both panels that those constraints restricts the X 3 
VEV u to u > 2 MeV. On the other hand, domain wall 
constraints [ H>j, EP] lead to the upper bound u < 10 
MeV, mentioned earlier. A noteworthy fact is that the 
allowed ranges of the mixing angles 0 12 , $13 — given in 
Table [T] only from relic density constraints — are fur¬ 
ther reduced to 4.5 x 10 -27 ^ $12 Ss 1.67 x 10 - 26 and 
1.0 x 10 ~ 13 ^ 0 i 3 ^ 2.75 x 10 -12 from the requirement 
of producing the correct X-ray and 7 -ray fluxes. The al¬ 
lowed ranges of the other parameters in Table [I] remain 
the same. 



u (MeV) u (MeV) 


Figure 8 : Allowed regions in the u — A 13 (left panel) and 
u — a (right panel) planes. 

In summary, our earlier model |36] can fit the anal¬ 
ysed data from Retll, while retaining the explanation 
for the 3.55 keV X-ray line — but with substantial mod¬ 
ifications. M x has to be pushed up to around 250 
GeV. Further, W + W~ need to replace bb among the 
decay products of x , 2 as the primary source of the 7 - 
excess. This new seed mechanism requires new Boltz¬ 
mann equations. They have been formulated with their 
consequences quantitatively worked out. The compati- 
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bility with other indirect constraints has been checked. 
The entire picture hangs together. 

* * * 
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